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A mathematically assisted reconstruction of the initial focus of the yellow 
fever outbreak in Buenos Aires (1871) 

M L Fernandez,^ M Otero,N Schweigmann,!^ H G SolarP" 


We discuss the historic mortality record corresponding to the initial focus of the yellow 
fever epidemic outbreak registered in Buenos Aires during the year 1871 as compared 
to simulations of a stochastic population dynamics model. This model incorporates the 
biology of the urban vector of yellow fever, the mosquito Aedes aegypti, the stages of 
the disease in the human being as well as the spatial extension of the epidemic outbreak. 
After introducing the historical context and the restrictions it puts on initial conditions and 
ecological parameters, we discuss the general features of the simulation and the dependence 
on initial conditions and available sites for breeding the vector. We discuss the sensitivity, 
to the free parameters, of statistical estimators such as: final death toll, day of the year 
when the outbreak reached half the total mortality and the normalized daily mortality, 
showing some striking regularities. The model is precise and accurate enough to discuss the 
truthfulness of the presently accepted historic discussions of the epidemic causes, showing 
that there are more likely scenarios for the historic facts. 


I. Introduction 

Yellow fever (YF) is a disease produced by an 
arthropod borne virus (arbovirus) of the family 
flaviviridae and genus Flavivirus. The arthropod 
vector can be one of several mosquitoes and the 
usual hosts are monkeys and/or people. Wild 
mosquitoes of genus Haemagogus, Sabetes and 
Aedes are responsible for the transmission of the 
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virus among wild monkeys, such as the Brown 
Howler Monkey (Alouatta guariba) associated to 
recent outbreaks of YF in Brazil, Paraguay and Ar¬ 
gentina [I]). In contrast, urban YF is transmitted 
by a domestic and anthropophilic mosquito, Aedes 
aegypti, human beings being the host [5]. Aedes ae¬ 
gypti is a tree hole mosquito, with origins in Africa, 
that has been dispersed through the world thanks 
to its association with people. 

During the end of the XVIII and the XIX cen¬ 
turies, YF caused large urban outbreaks in the 
Americas from Boston (1798), New York (1798) 
and Philadelphia (1793, 1797, 1798, 1799) in the 
North [3] to Montevideo (1857) and Buenos Aires 
(1858,1870, 1871) [1] in the South. These historical 
episodes arise as ideal cases for testing the capabili¬ 
ties of YF models in urban settings. Is it possible to 
reconstruct the evolution of one of these epidemic 
outbreaks? Can enough information be recovered 
to produce a thorough test on the models? This 
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is seldom the case, for example, for the study of 
the Memphis (1878) epidemic, with over 10000 ca¬ 
sualties, only 1965 were considered potentially us¬ 
able [S]. In contrast, the records of the outbreak in 
Buenos Aires 1871, unearthed and digitized for this 
work, left us with an amount of 1274 death cases 
located in time and space for the initial focus in 
the quarter of San Telmo, about 78% of the total 
mortality in the quarter [^. According to the 1869 
national census [7| Buenos Aires had 177787 inhabi¬ 
tants, 12329 of them living in San Telmo, about half 
of them just immigrated into the country mostly 
from Europe. 

In this work, we will compare the initial devel¬ 
opment of the epidemic outbreak (Buenos Aires, 
1871) with the simulations resulting from an eco- 
epidemiological model developed in Refs. [HHin], 
testing the worth of the predictive model. 

The simulations were performed under a number 
of assumptions, most of them essentially forced by 
the lack of better information. We will assume that: 

1. Now and before, YF is the same illness, i.e., 
we can use current information on YF devel¬ 
opment such as: the average extent of the in¬ 
cubation, infection, recovery, and toxic peri¬ 
ods, as well as the mortality level in 1871. In 
other words, the virus presents no substantial 
changes since 1871 to present days. We do 
not expect this hypothesis to be completely 
correct: the YF virus is an AAA-virus as op¬ 
posed to the stable £>AA-viruses, as such, mu¬ 
tations in about 140 years of continuous repli¬ 
cations in mosquitoes and primates can hardly 
be ruled out. Furthermore, present-day YF 
has been subject to different evolutionary pres¬ 
sures than the YF in the XIX century. While 
in the XIX century yellow fever circulated con¬ 
tinuously in human populations, today the 
wild part of the cycle involving wild popula¬ 
tions of monkeys plays a substantial role. 

2. The epidemic was transmitted by Aedes ae- 
gypti. There is no evidence of this fact since 
the scientific society and medical doctors in 
general were not aware of the role played by 
the mosquito until the confirmation given by 
Reed [3] of Finlay’s ideas [IT]. 0 

^According to other sources, it was Beauperthuy m 
the first one to accurately describe the transmission of YF 


We assume that Aedes aegypti has not changed 
since then, and/or there are no substantial 
changes in the life cycle, vector capabilities 
and adaptation between the (assumed) pop¬ 
ulation in 1871 and present-day populations 
in Buenos Aires city. After the eradication 
campaign (1958-1965) [13] . Aedes aegypti was 
eradicated from Buenos Aires m- Hence, the 
present populations result from a re-infestation 
and they are not the direct descendants of the 
mosquitoes population of 1871. 

3. Lacking time statistics for the duration of the 
different stages in the development of the ill¬ 
ness, reproduction of the virus and life cycle 
of the mosquito, we use, as distribution for 
such events, a maximum likelihood distribu¬ 
tion subject to the constrain of the average 
value for the cycle. In short, we use exponen¬ 
tially distributed times for the next event for 
all type of events. 

4. Finally, and most importantly, we assume that 
the human population mobility is not a factor 
in the local spread of the disease. We antic¬ 
ipate one of our conclusions: this assumption 
is likely to be false for the full development of 
the epidemic outbreak but seems reasonable 
for the early (silent) development. The study 
of the secondary foci of the epidemic outbreak 
merits a detailed analysis of the social and po¬ 
litical circumstances related to human mobility 
and it is beyond the possibilities of this study. 

Since we want the test to be as demanding as 
possible, more information is needed to simulate 
the outbreak eliminating sources of ambiguity and 
parameters to be fitted using the same test data. 
We recovered the following information: 

1. Estimations of daily temperatures. They are 
relevant since the temperature regulates the 
developmental rates of the mosquitoes. 

2. A very rough, anecdotal, estimation of the 
availability of breeding sites (BS) that, ulti¬ 
mately, control the carrying capacity of the 

by mosquitoes, as observed in the epidemic outbreak at 
Cumana, Venezuela (1853), as well as the efficient measures 
of protection taken by Native Americans, the use of nets to 
prevent the spread of the epidemic 1131 . 
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environment, the number of vectors and the 
infection rate. 

3. Human populations discriminated by block in 
the city. 

4. Estimations of the date of arrival of the virus 
to the city, putting bounds to the reasonable 
initial conditions for the simulation. 

This climatological, social and historical informa¬ 
tion represents a determining part of the recon¬ 
struction as it is integrated into the model jointly 
with the entomological and medical information 
to produce stochastic simulations of possible out¬ 
breaks to be compared with the historic records of 
casualties. 

We will show that the model predicts large prob¬ 
abilities for the occurrence of YF in the given his¬ 
torical circumstances and it is also able to answer 
why a minor outbreak in 1870 did not progress to¬ 
wards a large epidemic. The total number of deaths 
and the time-evolution of the death record will be 
shown to agree between the historical record and 
the simulated episodes as well, within the original 
focus. 

The rest of the manuscript will be organized as 
follows: we will begin with the description of YF in 
section II, including the eco-epidemiological model. 
In section III, we will address the relevant climato¬ 
logical, social and historical aspects. In section IV, 
we will explore the sensitivity of the model to ini¬ 
tial conditions and the number of available breed¬ 
ing sites, discussing the statistics more clearly influ¬ 
enced by vector abundance. The historic mortality 
records and the simulated records are compared in 
section |V]] We will finally discuss the performance 
of the model in section I VI. I 

II. The disease 

We will simply quote the fact sheet provided by the 
World Health Organization m as the standardized 
description: 

“YF is a viral disease, found in tropi¬ 
cal regions of Africa and the Americas. 

It principally affects humans and mon¬ 
keys, and is transmitted via the bite of 
Aedes mosquitoes. It can produce devas¬ 
tating outbreaks, which can be prevented 


and controlled by mass vaccination cam¬ 
paigns. 

The first symptoms of the disease usually 
appear 3-6 days after infection. The first, 
or acute, phase is characterized by fever, 
muscle pain, headache, shivers, loss of ap¬ 
petite, nausea and vomiting. After 3-4 
days, most patients improve and symp¬ 
toms disappear. However, in a few cases, 
the disease enters a toxicphase: fever 
reappears, and the patient develops jaun¬ 
dice and sometimes bleeding, with blood 
appearing in the vomit (the typical vomito 
negro). About 50% of patients who enter 
the toxic phase die within 10-14 days”. 

We add that the remission period lasts between 2 
and 48 hours [m, and as it was mentioned in the 
introduction, not only Aedes mosquitoes transmit 
the disease. 

i. The model 

The yellow fever model is rather similar to the al¬ 
ready presented dengue model m. the similarity 
corresponds to the fact that dengue is produced 
by a Flavivirus as well, it is transmitted by the 
same vector and follows the same clinical sequence 
in the human being, although with substantially 
lesser mortality. 

The model describes the life cycle of the mosquito 
[8] and its dispersalafter a blood meal, seeking 
oviposition sites [S]. The mosquito goes through 
several stages: egg, larva, pupa, adult (non parous), 
flyer (i.e., the mosquito dispersing) and adult 
(parous). In each stage, the mosquito can die or 
continue the cycle with a transition rate between 
the subpopulations that depends on the tempera¬ 
ture. The mortality in the larva stage is nonlinear 
and it regulates the population as a function of the 
availability of breeding sites. Thus, the transitions 
from adult to flyer are associated with blood meals, 
the event that can transmit the virus from human 
to mosquito and vice-versa. From the epidemiolog¬ 
ical point of view, the mosquito follows a SEI se¬ 
quence (Susceptible, Exposed —extrinsic period—, 
Infective). Correspondingly, the adult populations 
are subdivided according to their status with re¬ 
spect to the virus. We assume that there is no 
vertical transmission of the virus and eggs, larvae. 
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pupae and non parous adults are always suscepti¬ 
ble. 

The humans are subdivided in subpopulations 
according to their status with respect to the ill¬ 
ness as: susceptible (S), exposed (E), infective (I), 
in remission (r), toxic (T) and recovered (R). The 
temporary remission period is followed by recovery 
with a probability between 0.75 and 0.85 or a toxic 
period (probability 0.25 to 0.15) which ends half of 
the times in death and half of the times in recovery. 
The yellow fever model differs in the structure from 
the dengue model in Ref. m, as the human part 
of the dengue model is SEIR and the yellow fever 
model is SEIrRTD. However, the additional stages 
do not alter the evolution of the epidemic since the 
“in remission”, toxic and dead stages do not par¬ 
ticipate in the transmission of the virus. The YE 
parameters are presented in Table [T] 


Period 

value 

range 

Intrinsic Incubation (IIP) 

4 days 

3-6 days 

Extrinsic Incubation (EIP) 

10 days 

9—12 days 

Human Viremic (VP) 

4 days 

3-4 days 

Remission (rP) 

1 days 

0-2 days 

Toxic (tP) 

8 days 

7-10 days 


Probability 

value 

range 

Recovery after remission (rar) 

0.75 

0.75-0.85 

Mortality for toxic patients (mt) 0.5 


Transmission host to vector (ahv) 0.75 


Transmission vector to host (avh) 0.75 



Table 1: Parameters (mean value of state) adopted 
for YF. The range indicated is taken from PAHO 

na. 

The model is compartmental, all populations 
are counted as non-negative integers numbers and 
evolve by a stochastic process in which the time 
of the next event is exponentially distributed and 
the events compete with probabilities proportional 
to their rates in a process known as density- 
dependent-Poisson-process m- The model can be 
understood qualitatively with the scheme of the 
Fig. [H The model equations are summarized in 
Appendix 

The city is divided in blocks, roughly follow¬ 
ing the actual division (see Fig. [2D. The human 
populations are constrained to the block while the 
mosquitoes can disperse from block to block. 



Figure 1: Scheme of the yellow fever model. On the 
left side, the evolution of the mosquito and, on the 
right side, the evolution of human subpopulations. 
Hollow arrows indicate the progression through the 
life cycle of the mosquito following the sequence: 
egg, larva, pupa, adult (non-parous), flyer, adult 
(parous) and the repetition of the two last steps. 
The mortality events are not shown to lighten the 
scheme. Eggs are laid in the transition from flyer to 
adult. The adult mosquito populations are subdi¬ 
vided according to their status with respect to the 
virus as: susceptible (S), exposed (E) and infective 
(I). The virus is transmitted from mosquitoes to hu¬ 
mans and vice-versa in the transition from adult to 
flyer (blood meal) when either the mosquito or the 
human is infective and the other susceptible (red 
arrows). The red bold arrows indicate the progres¬ 
sion of the disease, from exposed to infective in the 
mosquito and, in humans, following the sequence: 
exposed (E), infective (I), in remission (r), toxic, 
recovered (R) or dead. 

III. Historical, social and climato¬ 
logical information 

i. When and how the epidemic started 

The YF outbreak in Buenos Aires (1871) was one 
of a series of large epidemic outbreaks associated 
to the end of the War of the Triple Alliance or 
Paraguayan war. The war confronted Argentina, 
Brazil and Uruguay (the three allies) on one side 
and Paraguay on the other side, and ended by 
March, 1870. By the end of 1870, Asuncion, 
Paraguay’s capital city, was under the rule of the 
Triple Alliance. The return of Paraguay’s war pris- 
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oners from Brazil (where YF was almost endemic 
at that time) to Asuncion triggered a large epi¬ 
demic outbreak [1] . The allied troops received their 
main logistic support from Corrientes (Argentina), 
a city with 11218 inhabitants according to the 1869 
census [7], located about 300 km south of Asun¬ 
cion (following the waterway) and 1000 km north 
of Buenos Aires along the Parana river (see Fig. 
131). On December 14, 1870, the first case of YF 
was diagnosed in Corrientes m. and a focus de¬ 
veloped around this case imported from Asuncion. 
According to some sources, the epidemic produced 
panic, resulting in about half the population leav¬ 
ing the city between December 15 and January 15 
m- However, other historical reasons might have 
played a relevant role since the city of Corrientes 
was under the influence and ruling of Buenos Aires, 
while in the farmlands, the General Ricardo Lopez 
Jordan was commanding a rebel army (a sequel of 
Argentina civil wars and the war of the Triple Al¬ 
liance). The subversion ended with the battle of 
Naembe, about 200 km east of Corrientes, on Jan¬ 
uary 26, 1871. 

Putting things in perspective, we must realize 
that in those times, YF was recognized only in its 
toxic stage associated to the black vomit, it is then 
perfectly plausible that recently infected individu¬ 
als would have left Corrientes and Asuncion reach¬ 
ing Buenos Aires, despite quarantine measures that 
were late and leaky o The death toll in Corri¬ 
entes was of 1289 people in the city (and about 700 
more in places around the city) |20| . representing 
a 11,5% of the population (notice that this num¬ 
ber is not consistent with current numbers in use 
by WHO [T7] which indicate a 7.5% of mortality 
in diagnosed cases of YF but is well in line with 
historical reports [23] of 20% to 70% mortality in 
diagnosed cases —the statistical basis has changed 
with the improved knowledge of early, not toxic, 
YF cases. 

According to this historical view, the initial ar¬ 
rival of infectious people to Buenos Aires happened, 
more likely, during December 1870 and January 
1871. In his study of the YF epidemic, written 
twenty three years after the epidemic outbreak, 
Jose Penna (MD) |H quotes the issue of the journal 

^On December 16, a sanitary official from Buenos Aires 
was commissioned to Corrientes to organize the quarantine, 
a measure that was applied to ships coming from Paraguay 
but not to those with Corrientes as departing port. 




Figure 2: Police districts from a map of the time 
and computer representation. The red colored area 
in district 14 (San Telmo) are the two blocks where 
the 1871 epidemic started. The green colored area 
in district 3 is the block where the Hotel Roma was 
and where the 1870 focus began (see section rVIlii.l) . 
The red and green lines sourround the region sim¬ 
ulated for the 1871 epidemics and the 1870 focus, 
respectively. Notice that districts 15 and 13 dis¬ 
agree in the maps. The computer representation 
follows the information in Ref. [21] from where the 
population information was obtained. 


“Revista Medico Quirtrrgica”, published in Buenos 
Aires on December 23, 1870 |21j, which presents a 
report regarding the sanitary situation during the 
last fifteen days, indicating the emergence of a “bil¬ 
ious fever” and a general tendency of other fevers 
to produce icterus or jaundice. In the next issue, 
dated January 8, 1871, the “Revista” indicates an 
important increase in the number of bilious fever 
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cases reported [35]. In a separate article, the doc¬ 
tors call the attention on how easily and how of¬ 
ten the quarantine to ships coming from Paraguay 
is avoided, and calls for strengthening the mea¬ 
sures. Penna indicates that the “bilious fever” (not 
a standard term in medicine) likely corresponded to 
milder cases of YF. We will term this idea “Penna’s 
conjecture” and will come back to it later. 

For our initial guess, we considered this informa¬ 
tion as evidence that the epidemic outbreak started 
during December, 1870. Exploring the model, and 
arbitrarily, we took December 16, 1870, as the 
time to introduce two infectious people with YF in 
the simulations, at the blocks where the mortality 
started. Yet, an educated guess for Penna’s con¬ 
jecture is to consider the 3-6 days needed from in¬ 
fection to clinical manifestation and the 9-12 days 
of the extrinsic cycle. Hence, since the first clini¬ 
cal manifestations of transmitted YF happened be¬ 
tween December 11-23, we would guess the infected 
people arriving somewhere between November 21 
and December 11. 


Bolivia ^ 



Figure 3: A 1870 map |2S] showing Asuncion next 
to the Paraguay river, Corrientes and Rosario next 
to the Parana river and Buenos Aires (spelled 
Buenos Ayres) next to the Rio de la Plata. 

Yet, we must take into account that Penna’s con¬ 
jecture contrasts with the conjectures presented by 
MD Wilde and MD Mallo, members of the San¬ 


ity committee in charge during the YF epidemic. 
Wilde and Mallo advocated for the spontaneous 
origin of the disease, very much in line with the 
theories of miasmas in use in those times, theo¬ 
ries that guided the sanitary measures taken m- 
Wilde and Mallo also argued that Asuncion could 
not be the origin of the epidemic, because of their 
belief that the ten or fifteen days quarantine (count¬ 
ing since the last port touched) was enough to avoid 
the propagation of the disease. This belief contrasts 
with the experience of 1870 (in Buenos Aires) where 
a ten day quarantine was not enough to prevent a 
minor epidemic [3]. Nevertheless, the quarantine 
measures werefully implemented in Corrientes by 
December 31, 1870. The measures were later lifted 
because of the epidemic in Corrientes and imple¬ 
mented at ports down the Parana river, being com¬ 
pleted nearby Buenos Aires (ports of La Conchas, 
Tigre, San Fernando and “La Boca” within Buenos 
Aires city) by mid-February when the epidemic was 
in full development in Buenos Aires according to 
the port sanitary authorities, Wilde and Mallo m- 

Being Corrientes the source of infected people 
can hardly be disregarded. With about 5000 people 
leaving the city between December 15 and January 
15 [T^, a city where YF was developing. 

According to Wilde and Mallo nni, there were 
(non-fatal?) YF cases in Buenos Aires as early 
as January 6, 1871, reported by MDs Argerich 
and Gallarani as well as documented cases of YF 
death after disembarking in Rosario (200 km North 
of Buenos Aires along the Parana river) having 
boarded in Corrientes. 

ii. Breeding sites 

One of the key elements in the reconstruction 
and simulation of an epidemic transmitted by 
mosquitoes is to have an estimation of their num¬ 
bers which will be reflected directly in the propa¬ 
gation of the epidemic. In the mosquito model |H| , 
this number is regulated by the quality and abun¬ 
dance of breeding sites. The production of a single 
breeding site, normalized to be a flower pot in a 
local cemetery, was taken as unit in Ref. |S] and 
the number of breeding sites measured in this unit 
roughly corresponds to half a liter of water. 

The Aedes aegypti population in monitored areas 
of Buenos Aires, today, is compatible with about 
20 to 30 breeding sites per block [Hj. Estimating 
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the number of sites available for breeding today is 
already a difficult task, the estimation of breeding 
sites available in 1871 is a nearly impossible one. 
In what remains of this subsection, we will try to 
get a very rough a-priori estimate. 


District 

Population/(100 m)^ 

BS/(100 m)-^ 

1 

339 

391.0 

2 

279 

300.0 

3 

428 

522.0 

4 

353 

443.0 

5 

330 

430.0 

6 

259 

365.0 

13 

160 

196.0 

14 

224 

300.0 

15 

90 

157.0 

16 

165 

316.0 

18 

23 

52.0 

19 

13 

26.0 

20 

30 

52.0 


Table 2: Population data. Buenos Aires, 1869 [7]. 
Population density by police district (see Fig. ID) 
and equivalent breeding sites, BS, originally esti¬ 
mated as proportional to the house density in the 
police district. 

A very important difference between those days 
and the present corresponds to the supply of fresh 
water which today is taken from the river, pro¬ 
cessed and distributed through pipes; but in those 
days, it was an expensive commodity taken from 
the river by the “waterman” and sold to the cus¬ 
tomers who, in turn, had to let it rest so that the 
clay in suspension decanted to the bottom of the 
vessel (a process that takes at least 3 days). Ad¬ 
ditionally, there were some wells available but the 
water was (is) of low quality (salty). The last, and 
rather common resource mm, was the collection 
of rain water in cisterns. 

iii. Temperature reconstruction 

Aedes aegypti developmental times depend on tem¬ 
perature. Although it would seem reasonable to 
use as substitute of the real data of the average 
temperature registered since systematic data collec¬ 
tion began, records of temperature in those times 
were kept privately m and are available. The data 
set consists of three daily measurements made from 
January 1866 until December 1871, at 7AM, 2PM 
and 9PM. When averaged, the records allow an es¬ 
timation of the average temperature of the day bet- 


Januaiy 1st, 1866 January 1st, 1871 



Time / days 


Figure 4: Average daily temperature and periodic 
approximation fitted accordng to Eq. o, t = 0 
corresponds to January 1, 1866. The inset shows 
the difference between the measured temperatures 
and the fit (residuals) during 1871. 

ter than the usual procedure of adding maximum 
and minimum dividing by two. Unfortunately, the 
register has some important missing points during 
the epidemic outbreak. Because of this problem, 
the data in Ref. m was used to fit an approxima¬ 
tion in the form: 

T =7.22°C X cos(27rt/(365.25 days) -f 5.9484) 

-h 17.21°C, (1) 

following Ref. [28] and then extrapolating to the 
epidemic period. In Fig. [U the data and the fit are 
displayed. The residuals of the fit do not present 
seasonality or sistematic deviations, as we can see 
in the inset of Fig. 

It is worth noticing that a similar fit on temper¬ 
ature data from the period 1980-1990 presents a 
mean temperature of 18.0°C, amplitude of 6.7°C 
and a phase shift of 6.058°C |S| (notice that t = 0 
in the reference corresponds to July 1 while in this 
work it corresponds to January 1). According to 
the threshold computations in Ref. |F|, the cli¬ 
matic situation was less favorable for the mosquito 
in 1866-1871 than in the 1960-1991 period. 

The reconstruction of temperatures needs to be 
performed at least from the 1868 winter, since a 
relatively arbitrary initial condition in the form of 
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eggs for July 1, 1868 is used to initialize the code, 
and then run over a transitory of two years. Such 
procedure has been found to give reliable results 
[5]. There are several factors in the biology of Ae. 
ae. that indicate that the biological response to 
air-temperature fluctuations is reflected in attenu¬ 
ated fluctuations of biological variables. First, the 
larvae and pupae develop in water containers, thus, 
what matters is the water temperature. This fact 
represents a first smoothing of air-temperature fluc¬ 
tuations. Second, insects developmental rates for 
fluctuating temperature environments correspond 
to averages in time of rates obtained in constant 
temperature environments m, an alternative view 
is that development depends on accumulated heat 
[30]. Such averages occur over a period of about 6 
days at 30 °C and longer times for other tempera¬ 
tures and non-optimal food conditions |31j . Third, 
the biting rate (completion of the gonotrophic cy¬ 
cle) depends as well on temperatures averaged over 
a period of a few days. Last, mosquitoes actively 
seek the conditions that fit best to them and more 
often than not, they are found resting inside the 
houses. 

iv. Mortality data 

For this work, the daily mortality data recorded 
during the 1871 epidemic outbreak |5] is key. This 
statistical work has received no attention in the 
past, and no study of the YF outbreak in Buenos 
Aires made reference to this information. We have 
cross-checked the information with data in the 1869 
national census |7], as well as with data in published 
works |3|, and the details are consistent among 
these sources. 

The data set is presented here closing the historic 
research part. 

V. Revising the clinical development of yel¬ 
low fever. 

A YF epidemic outbreak happened in Buenos 
Aires, in 1870, developing about 200 cases ji] (the 
text is ambiguous on whether the cases are toxic 
or fatal). The epidemic outbreak was noticed by 
February 22 (first death), a sailor who left Rio de 
Janeiro (Brazil) on February 7, and presumably 
landed on February 17 (no cases of YF were re¬ 
ported on board of the Poitou —the boat). 



Figure 5: Daily mortality cases in the police district 
14 (see Fig. jj) corresponding to the quarter San 
Tehno |5], t = 0 corresponds to January 1, 1871. 

This well documented case allows us to see the 
margins of tolerance that have to be exercised in 
taking medical information prepared for clinical 
use as statistical information. Assume, following 
Penna, that the sailor was exposed to YF before 
boarding in Rio de Janeiro, according to informa¬ 
tion in Table |T] collected from the Pan American 
Health Organization El, adding incubation and 
viremic period, we have a range of 6-10 days, hence 
the sailor was close to the limit of his infectious pe¬ 
riod. He was not toxic, according to the MD on 
board who signed a certihcate accepted by the san¬ 
itary authority. Yet, five days later, he was dying, 
making the remission plus toxic period of 5 days, 
shorter than the range of 7-12 days listed in Ta- 
bledjand substantially shorter than the 10-15 days 
(remission plus toxic) communicated in Ref. |16j . 
Shall we assume as precise the values reported in 
Ref. [161[T7] we would have to conclude that the 
disease was substantially different, at least in its 
clinical evolution, in 1871 as compared with present 
days. 

In clinical studies performed during a YF epi¬ 
demic on the Jos Plateau, Nigeria, Jones and Wil¬ 
son [221 report a 45.6% overall mortality and a 
significant difference in the duration of the illness 
for fatal and non-fatal cases with averages of 6.4 
and 17.8 days. Serie et al. [SSj reports for the 
1960-1962 epidemic in Ethiopia a mortality rang- 
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ing from 43% (Koure) to 100% (Boloso) and 50% 
(Menera) with a total duration of the clinic phase 
of the illness of 7.14, 2.14 and 4.5 days, respectively 
(weighted average of 4.6 days in 18 cases). 

We must conclude that the extension of the toxic 
period preceding the death presents high variabil¬ 
ity. This variability may represent variability in the 
illness or in medical criteria. For example, Serie 
[55] indicates that the 100% mortality found at the 
Boloso Hospital is associated to the admission cri¬ 
teria giving priority to the most severe cases. In 
correspondence with this extremely high mortality 
level, the survival period is the shortest registered. 
The minimum length of the clinical phase is of 10- 
14 days or 13-18 days, depending of the source 
(adding viremic, remission and toxic periods). We 
note that not only the toxic period of fatal cases 
must be shorter than the same period for non-fatal 
cases, but also the viremic period must be shorter 
in average, if all the pieces of data are consistent. 

The time elapsed between the first symptoms and 
death is probably longer today than in 1871, since 
it, in part, reflects the evolution of medical knowl¬ 
edge. The hospitalization time is also rather arbi¬ 
trary and changes with medical practices which do 
not reflect changes in the disease. 

A rudimentary procedure to correct for this dif¬ 
ferences is to shift the simulated mortality some 
fixed time between 5 and 8 days (the difference be¬ 
tween our 13 days guessed (Table [T]) and the 4.5- 
6.4 reported for Africa 151155I). Such a procedure 
is not conceptually optimal, but it is as much as 
it can be done within present knowledge. We cer¬ 
tainly do not know whether just the toxic period 
must be shortened or the viremic period must be 
shortened as well, and in the latter case, how this 
would affect the spreading of the disease. 

A second source of discrepancies between 
recorded data and simulations are the inaccuracies 
in the historic record. Can we consider the daily 
mortality record as a perfect account? Which was 
the protocol used to produce it? We can hardly 
expect it to be perfect, although we will not make 
any provision for this potential source of error. 

IV. Simulation results 

The simulations were performed using a one-block 
spatial resolution, with the division in square blocks 


of the police districts 14 (San Telmo), 16, 2 and 
4, corresponding to Concepcion, Catedral Sur and 
Montserrat; and part of the districts 6, 18, 19, 19A 
and 20 (see Fig. H]), and totalized for each po¬ 
lice district to obtain daily mortality comparable 
to those reported in Ref. (5] and picture in Fig. 

Numerical mosquitoes were not allowed to fly 
over the river. At the remaining borders of the 
simulated region, a Stochastic Newmann Bound¬ 
ary Condition was used, meaning that the mosquito 
population of the next block across the boundary 
was considered equal to the block inside the region; 
but the number of mosquito dispersion events as¬ 
sociated to the outside block was drawn randomly, 
independently of the events in the corresponding in¬ 
ner block. Larger regions for the simulations were 
tested producing no visible differences. 

The time step was set to the small value of 30 
s, avoiding the introduction of further complica¬ 
tions in the program related to fast event rates for 
tiny populations ED, although an implementation 
of the method in Ref. |51j , not relying on the small¬ 
ness of the time step so heavily, is desirable for a 
production phase of the program. 

Before we proceed to the comparison between 
the historic mortality records of the epidemic and 
the simulated results, we need to gain some un¬ 
derstanding regarding the sensitivity of the simula¬ 
tions to the parameters guessed and the best forms 
of presenting these results. We performed a mod¬ 
erate set of computations, since the code has not 
been optimized for speed and it is highly demand¬ 
ing for the personal computers where it runs for 
several days. Here, we illustrate the main lessons 
learned in our explorations. Poor people, unable to 
buy large quantities of water, had to rely mostly 
on the cisterns and other forms of keeping rain wa¬ 
ter. Since 1852, when the population of Buenos 
Aires was about 76000 people, there was an impor¬ 
tant immigration flow, increasing the population 
to about 178000 people by 1869 l^. The immi¬ 
grants occupied large houses where they rented a 
room, usually for an entire family, a housing that 
was known as “conventillo” and was the dominant 
form of housing in some districts such as San Telmo, 
where the epidemic started [50]. In some police 
chronicles of the time, houses with as many as 300 
residents are mentioned m- Under such difficult 
social circumstances, we can only imagine that the 
number of breeding sites available to mosquitoes 
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has to be counted as orders of magnitude larger 
than present-day available sites. 

An a-priori and conservative estimation is to con¬ 
sider about ten times the number of breeding sites 
estimated today. Thus, we assume, as a first guess, 
300 breeding sites per block in San Telmo. We 
will have to tune this number later as it regulates 
mosquito populations and the development of the 
epidemic focus. The number of breeding sites is the 
only parameter tuned to the results in this work. 

More precisely, the criteria adopted was to make 
the number of (normalized) BS proportional to the 
number of houses per block taken from historic 
records EH, adjusting the proportionality factor to 
the observed dynamics. We introduce the notation 
BSxY to indicate a multiplicative factor of Y. 

The population of each police district was set 
to the density values reported in the 1869 census 
and the police districts geography was taken 
from police records m and referenced according 
to maps of the city at the time I23S0]- 

Table E] shows the average population per block, 
initially estimated number of breeding sites and 
number of houses for the district of the initial focus, 
San Telmo ^/Id and nearby districts ( 7 / 16, 4, 2). 
A sketch of Buenos Aires police districts according 
to a 1887 map [37j is displayed alongside with the 
computer representation in Fig. [5] 

It is a known feature of stochastic epidemic mod¬ 
els m that the distribution of totals of infected 
people has two main contributions. One is that the 
small epidemic outbreaks when none or a few sec¬ 
ondary cases are produced and the extinction time 
of the outbreak comes quickly. The otheris that the 
large epidemic outbreaks which, if the basic repro¬ 
ductive number is large enough, present a Gaussian 
shape separated by a valley of improbable epidemic 
sizes from the small outbreaks. 

While the present model does not fall within the 
class of models discussed in Ref. IS], the gen¬ 
eral considerations applied to stochastic SIR mod¬ 
els qualitatively apply to the present study. Yet, 
simulations started early during the summer season 
follow the pattern just described in Ref. [Sj , but 
simulations started later do not present the proba¬ 
bility valley between large and small epidemics. 

We have found useful to present the results dis¬ 
aggregated in the form: epidemic size, daily per¬ 
centage of mortality relative to the total mortality 
and time to achieve half of the final mortality. This 


presentation will let us realize that most of the fluc¬ 
tuation is concentrated in the total epidemic size, 
while the daily evolution is relatively regular, ex¬ 
cept, perhaps, in the time taken to develop up to 
50% of the mortality (depending on the abundance 
of vectors and the initial number of infected hu¬ 
mans and chance). 

i. Total mortality (epidemic size) 

Since historical records include mostly the number 
of causalities, it appears sensible for the purposes 
of this study to use the total number of deaths as 
a proxy statistics for epidemic size. 



Total mortality 


Figure 6 : San Telmo. Total mortality histograms 
for different number of breeding sites, computed af¬ 
ter 100 simulations with the same initial condition 
corresponding to 2 infectious people located in San 
Telmo on January 1, 1871, at the same location 
where the initial death happened in the historical 
event. From top to bottom, multiplication factors 
(bin-width : frequency of no-epidemic) BSxl (126.2 
: 0.23), BSx2 (209.8 : 0.09), BSx3 (128.2 : 0.06) 
and BSx4 (50.6 : 0.02). The y-axis indicates fre¬ 
quency in a set of 100 simulations. 

The total mortality depends strongly on the 
stochastic nature of the simulations, initial condi¬ 
tions and ecological parameters guessed. Qualita¬ 
tively, the results agree with the intuition, although 
this is an a-posteriori statement, i.e., only after see¬ 
ing the results we can find intuitive interpretations 
for them. 
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The discussion assumes that the development of 
the epidemic outbreak was regulated by either the 
availability of vectors (mosquitoes) or the exhaus¬ 
tion of susceptible people, the first situation repre¬ 
sents a striking difference with standard SIR mod¬ 
els without seasonal dependence of the biological 
parameters. Actually, in Fig. we compare fre¬ 
quencies of epidemics binned in five bins by final 
epidemic size for different sets of 100 simulations 
with different number of breeding sites. The num¬ 
ber of breeding sites is varied in the same form 
all along the city, keeping the proportionality with 
housing, and it is expressed as multiplicative factor 
(BSxl=l, BSx2=2, BSx3=3, BSx4=4) presented 
in Table [5] Notice also that the width of the bins 
progresses as 126.2, 209.8, 128.2 and 50.6 indicat¬ 
ing how the dispersion of final epidemic sizes first 
increases with the number of breeding sites but for 
larger numbers decreases. 

We can see how for a factor 2 (BSx2) (and larger) 
the mortality saturates, indicating the epidemic 
outbreaks are limited by the number of susceptible 
people. For our original guess, factor 1, the epi¬ 
demic is limited by the seasonal presence/absence 
of vectors. However, for higher factors, there is 
a substantial increase in large epidemic outbreaks 
with larger probabilities for larger epidemics. Only 
for the factor 4 (BSx4) the most likely bin includes 
the historical value of 1274 deaths. 

A second feature, already shown in the Dengue 
model m, is that outbreaks starting with the ar¬ 
rival of infectious people in late spring will have a 
lesser chance to evolve into a major epidemic. Yet, 
those that by chance develop are likely to become 
large epidemic outbreaks since they have more time 
to evolve. On the contrary, outbreaks started in 
autumn will have low chances to evolve and not 
a large number of casualties. The corresponding 
histograms can be seen in Fig. [71 Figure [7] also 
shows how the outbreaks that begin on December 
16, as well as simulations starting on January 1, 
present higher probabilities of large epidemics than 
of small epidemics, but this tendency is reverted 
in simulations of outbreaks that start by February 
16. This transition is, again, the transition between 
outbreaks regulated by the number of available sus¬ 
ceptible humans and those regulated by the pres¬ 
ence or absence of vectors. 



Figure 7: San Telmo. Total mortality histograms 
for different date of arrival of infectious people. 
Computed after 100 simulations with the same ini¬ 
tial condition corresponding to 2 infectious people 
located in San Telmo and a density factor BSx2.5. 
From top to bottom: December 16, 1870; January 
1, January 16 and February 15, 1871. 

ii. Mortality progression 

One of the most remarkable facts unveiled by the 
simulations is that when the time evolution of the 
mortality is studied as a fraction of the total mor¬ 
tality, much of the stochastic fluctuations are elimi¬ 
nated and the curves present only small differences 
(see Fig. [8]). The similarity of the normalized evo¬ 
lution allows us to focus on the time taken to pro¬ 
duce half of the mortality (labelled Txji)- 

We notice, in Fig. |8l that the T 1/2 in these runs 
lay between 69-110, compared to the historic values 
ofTi/2 =73. 

The latter observation brings the attention to a 
remarkable fact of the simulations: not only the 
normalized progression of the outbreaks are rather 
similar but also, there is a correspondence between 
early development and large mortality. Drawing 
Ti /2 against total-mortality (Fig. [9]), we see that 
even for different number of breeding sites, all the 
simulations indicate that the final size is a noisy 
function of the day when the mortality reaches half 
its final value. This function is almost constant for 
small Ti /2 and becomes linear with increasing dis¬ 
persion when T 2/2 is relatively large. Once again, 
the two different forms the outbreak is controlled. 


050002-11 



















Papers in Physics, vol. 5, art. 050002 (2013) /ML Fernandez et al. 



— average 

— average - variance 

— average + variance 


Day (1 is January 1st 1871) 


Figure 8: Evolution of the mortality level for all the 
epidemics with three or more fatal cases in a batch 
of 100 simulations (96 runs) differing only in the 
pseudo-random series. The breeding sites number 
has a factor 4, and the initial contagious people 
were placed on January 1, 1871, in San Telmo. 
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Figure 9: San Telmo. Total mortality versus day 
when the mortality reached half its final number for 
different number of breeding sites. The mortality 
appears to be mainly a function of T 1/2 and roughly 
independent of the number of breeding sites. 

V. Real against simulated epidemic 

We would like to establish the credibility of the 
statement: the historic mortality record for the San 
Telmo focus belongs to the statistics generated by 


the simulations. 

To achieve this goal, we need to compare the 
daily mortality in the historical record and the sim¬ 
ulations. Since, in the model, the mortality pro¬ 
ceeds day after day with independent random in¬ 
crements (as a consequence of the Poisson charac¬ 
ter of the model), it is reasonable to consider the 
statistics 


2 _ f HM{i)- MM{i) 

where i runs over the days of the year, HM{i) is 
the fraction of the total death toll in the historic 
record for the day i, MM{i) is the average of the 
same fraction obtained in the simulations and D{i) 
is the corresponding standard deviation for the sim¬ 
ulations. The sum runs over the number of days in 
which the variance is not zero, for BSx3 and BSx4 
in no case D{i) = 0 and {HM{i) — MM(i) ^ 0). 
The number of degrees corresponds to the number 
of days with non-zero mortality in the simulations 
minus one. The degree discounted accounts for the 
fact that HM{i) = MM{i) = 1. 

i. Tuning of the simulations 

Before we proceed, we have to find an acceptable 
number of breeding sites, a reasonable day for the 
arrival of infected individuals (assumed to be 2 in¬ 
dividuals arbitrarily) and adjust for the uncertainty 
in survival time. Actually, moving the day of arrival 
d days earlier and shortening the survival time by 
d will have essentially the same effect on the sim¬ 
ulated mortality (providing d is small), which is to 
shift the full series by d days. This is, assigning to 
the day i the simulated mortality SM oi the day 
i + d, {SMii + d)). Hence, only two of the param¬ 
eters will be obtained from this data. 

As we have previously observed, the total mor¬ 
tality presents a large variance in the simulations. 
Moreover, in medical accounts of modern time 
[^133) . the mortality ranges between 46% and 
100 % while in historical accounts the percentage 
goes from 20 to 70 [53]. Hence, a simple adjustment 
of the mortality coefficient from our arbitrary 50% 
within such a wide range would suffice to eliminate 
the contributions of the total epidemic size. The 
average simulated epidemic for BSx4 is of ~ 1248 
deaths while the historic record is of 1274 deaths. 
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BSxY 

d 


degree 

Probability 

2 

8 

490.2 

177 

0 . 

3 

7 

191.7 

173 

0.15 

4 

3 

151.9 

168 

0.80 

4 

4 

143.8 

168 

0.91 

4 

5 

145.6 

168 

0.89 


Table 3: calculations according to ©■ We indi¬ 

cate the multiplicative factor applied to the breed¬ 
ing sites described in Tabled the shift applied to 
the statistics produce with the parameters of Table 
[U the value of the statistic , the number of de¬ 
grees of freedom, and the probability P{x > x^) for 
a random variable x distributed as a x^-distribution 
with the indicated degrees of freedom. The stan¬ 
dard reading of the x^ test indicates that, for BSx3 
and BSx4, the statement “the deviations from the 
simulations mean of the historic record (deviations 
assumed to be distributed as x^ with the indicated 
degrees) belong to the simulated set” is not found 
likely to be false by the test in the cases BSx3 and 
BSx4. 

Hence, to match the mean with the historic record, 
it would suffice to correct the mortality from 50% 
to 51%. 

We can disregard the idea that the epidemic 
started before December 20, 1870 (Penna’s con¬ 
jecture), since it is not possible to simultaneously 
obtain an acceptable final mortality and an accept¬ 
able evolution of the outbreak. Our best attempt 
corresponds to an epidemic starting by December 
14, 1870, which averages rs 1212 deaths and with a 
deviation of the mortality of x^ = 219.1 with 175 
degrees, giving a probability P{x > 219,1) = 0.01. 

We focus on arrival dates around January 1, 
1871. We can also disregard, for this initial condi¬ 
tion, the original guess BSxl corresponding to 300 
breeding sites per block in San Telmo, since it pro¬ 
duces too small epidemics. We present results for 
the epidemics corresponding to BSx2, BSx3 and 
BSx4 in Table [3] for different numbers of BS and 
shift d. 

ii. Comparison 

The conclusion of the x^ tests is that the simula¬ 
tions performed with BSx3 and beginning between 
December 24 1870 and January 1, 1871 (with a sur¬ 


vival period shortened between 0 and 8 days) are 
compatible with the historical record. However, the 
compatibility is larger when BSx4 is considered and 
the beginning of the epidemic is situated between 
December 28, 1870 and January 5, 1871 (with a 
survival period shortened between 0 and 8 days re¬ 
spectively). We illustrate this comparison with Fig. 

uni 



Days of 1871 

Figure 10: The historic record of accumulated mor- 
tality as a function of the day of the year 1871 is 
presented. The averaged accumulated mortality as 
well as curves shifted one standard deviation are 
shown for comparison. An acceptable look alike in¬ 
dividual simulation chose by visual inspection from 
a set of 100 simulations is included as well. Two 
exposed individuals (not yet contagious) were in¬ 
troduced in the same blocks where the historic epi¬ 
demic started the day January 1, 1871, with BSx4. 
The simulated mortality is anticipated in 4 days. 
Statistical estimates were taken as averages over 
96 runs which resulted in secondary mortality, out 
of a set of 100 runs. 


iii. The 1870 outbreak 

The records for the 1870 outbreak are scarce. Of 
the recognized cases, only 32 entered the Lazareto 
(hospital) and 19 of them were originated in the 
same block that the first case. Secondary cases are 
registered at the Lazareto’ books starting on March 
30 (two cases) and continuing with daily cases, the 
final outcome of these cases and the cases not en¬ 
tered at the Lazareto [3] are not clear. 
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Except for the precise initial condition, corre¬ 
sponding to one vireinic (infective) person located 
at the Hotel Roma (district 4 in Fig. [5]), the in¬ 
formation is too imprecise to produce a demanding 
test for the model. 

We performed a set of 100 simulations introduc¬ 
ing one viremic (infective) person at the precise 
block where the Hotel Roma was, on February 17. 
The number of breeding sites was kept at the same 
factor 4 with respect to the values tabulated in Ta¬ 
ble [5] that was used for the best results in the study 
of the 1871 outbreak. Needless to say, this does not 
need to be true, as the number of breeding sites may 
change from season to season. 

The distribution of the final mortality is shown 
in Fig. [TT] As we can see, relatively small epi¬ 
demics of less than 200 deaths cannot be ruled out, 
although there are much larger epidemics also likely 
to happen. 



Total mortality 


Figure 11: Mortality distribution for the simula¬ 
tions beginning on February 17 with the incorpo¬ 
ration of one infective (viremic) in the block of the 
Hotel Roma and BSx4. The histogram is the result 
of 87 runs which resulted in epidemics (13 runs did 
not result in epidemics). The width of the bins is 
322.2, and the first epidemic bin goes from 17 to 
339 deaths with a frequency of 15/100. 

Actually, a slow start of the epidemic outbreak 
would favor a small final mortality, as it can be seen 
in Fig. [121 Not only there is a relation between a 
low mortality early during the outbreak (such as 
April 15) and the final mortality, but we also see 


that the sharp division between small and large epi¬ 
demics is not present in this family of epidemic out¬ 
breaks differing only in the pseudo-random number 
sequence. 



Early mortality 


Figure 12: Final mortality against early mortality 
for two dates: March 30 and April 15. A low early 
mortality “predicts” a low final mortality as the out¬ 
break does not have enough time to develop. Sim¬ 
ulations correspond to the conditions of the 1870 
small epidemic. BSx4, one infected arriving to the 
Hotel Roma on February 17. The historic informa¬ 
tion indicates that secondary cases were recorded 
by March 30. Hence, corresponding to a slow start. 
The final mortality is not known, but it is believed 
it has been in the 100-200 range. 

For the smallest epidemics simulated, the sec¬ 
ondary mortality starts after March 30. Hence, the 
1870 focus can be understood as a case of relative 
good luck and a late start within more or less the 
same conditions than the outbreak of 1871. 

VI. Conclusions and final discussion 

In this work, we have studied the development of 
the initial focus in the YF epidemic that devastated 
Buenos Aires (Argentina) in 1871 using methods 
that belong to complex systems epistemology [12] . 

The core of the research performed has been the 
development of a model (theory) for an epidemic 
outbreak spread only by the mosquito Aedes ae- 
gypti represented according to current biological 
literature such as Christophers |2] and others. The 
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translation of the mosquito’s biology into a com¬ 
puter code has been performed earlier Pl^l and the 
basis for the spreading of a disease by this vector 
has been elaborated in the case of Dengue previ¬ 
ously m- The present model for YF is then an 
adaptation of the Dengue model to the particular¬ 
ities of YF, and the present attempt of validation 
(failed falsification) reflects also on the validity of 
this earlier work. 

The present study owes its existence to the work 
of anonymous police officers [6] that gathered and 
recorded epidemics statistics during the epidemic 
outbreak, in a city that was not only devastated by 
the epidemic, but where the political authorities 
left in the middle of the drama as well m- 

We have gathered (and implemented in a model) 
entomological, ecological and medical information, 
as well as geographic, climatological and social in¬ 
formation. After establishing the historical con¬ 
straints restricting our attempts to simulate the 
historical event, we have adjusted the density of 
breeding sites to be the equivalent to 1200 half¬ 
liter pots as those encountered in today Buenos 
Aires cemeteries (the number corresponds to San 
Telmo quarter). Perhaps a better idea of the num¬ 
ber of mosquitoes present is given by the maxi¬ 
mum of the average number of bites per person 
per day estimated by the model, which results in 5 
bites/(person day) (to be precise: the ratio between 
the maximum number of bites in a block during a 
week and the population of the block, divided by 
seven). 

The population of the domestic mosquito Aedes 
aegypti in Buenos Aires 1870-1871 was large 
enough to almost assure the propagation of YF 
during the summer season. The only effective 
measures preventing the epidemic were the natu¬ 
ral quarantine resulting from the distance to trop¬ 
ical cities were YF was endemic (such as Rio de 
Janeiro) and the relatively small window for large 
epidemics, since the extinction of the adult form 
of the mosquito during the winter months prevents 
the overwintering of YF. In this sense, the relatively 
small outbreak of 1870 is an example of how a late 
arrival of the infected individual combined with a 
touch of luck produced only a minor sanitary catas¬ 
trophe. 

By 1871, as a consequence of the end of the 
Paraguayan war and the emergence of YF in Asun¬ 
cion, the conditions for an almost unavoidable epi¬ 


demic in Buenos Aires were given. The interme¬ 
diate step taken in Corrientes, with the panic and 
partial evacuation of the city, adding the lack of 
quarantine measures, was more than enough to 
make certain the epidemic in Buenos Aires. On the 
contrary, Penna’s conjecture of an earlier starting 
during December, 1870 are inconsistent with the 
biological and medical times as implemented in the 
model. We can disregard this conjecture as highly 
improbable. 

The historic mortality record is consistent with 
an epidemic starting between December 28 and 
January 5, being the symptomatic period (viremic 
plus remission plus toxic) of the illness between 13 
and 5 days. Furthermore, the existence of iron- 
fatal cases of YF by January 6 mentioned by some 
sources m would be consistent, provided the cases 
were imported. 

In retrospect, the present research began as an 
attempt to validate/falsificate the YF model and, 
in more general terms, the model for the transmis¬ 
sion of viral diseases by Aedes aegypti using the 
historic data of this large YF epidemic. As the re¬ 
search progressed, it became increasingly evident 
that the model was robust. In successive attempts, 
every time the model failed to produce a reasonable 
result, it forced us to revise the epidemiological and 
historical hypotheses. In these revisions, we ended 
up realizing that the accepted origin of the epidemic 
in imported cases from Brazil, actually hides the 
central role that the epidemics in Corrientes had, 
and the gruesome failure of not quarantining Cor¬ 
rientes once the mortality started by December 16, 
1870, about two weeks before the deducted begin¬ 
ning of the outbreak in Buenos Aires. 

The same study of inconsistencies between the 
data and the reconstruction made us focus on the 
survival time of those clinically diagnosed with YF 
that finally die. The form in which the illness 
evolves anticipates the final result. Jones and Wil¬ 
son m indicate the symptoms of cases with a bad 
prognosis including the rapidity and degree of jaun¬ 
dice. This information suggests, in terms of model¬ 
ing, that death is not one of two possible outcomes 
at the end of the “toxic period”, as we have first 
thought. Separation of to-recover and to-die sub¬ 
populations could (should?) be performed earlier 
in the development of the illness, each subpopula¬ 
tion having its own parameters for the illness. Yet, 
while in theory this would be desirable, in prac- 
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tice it would have, for the time being, no effect, 
since the characteristic periods of the illness have 
not been measured in these terms. 

Epidemics transmitted by vectors come to an end 
either when the susceptible population has been 
sufficiently exposed so that the replication of the 
virus is slowed down (the classical consideration in 
SIR models) or when the vector’s population is dec¬ 
imated by other (for example, climatic) reasons. 
The model shows that both situations can be dis¬ 
tinguished in terms of the mortality statistics. 

We have also shown that the total mortality of 
the epidemic is not difficult to adjust by changing 
the death probability of the toxic phase, and as 
such, it is not a demanding test for a model. The 
daily mortality, when normalized, shows sensitivity 
to the mosquito abundance, specially in the evo¬ 
lution times involved, since the general qualitative 
shape appears to be fixed. In particular, the date in 
which the epidemic reaches half the total mortality 
is advanced by larger mosquito populations. How¬ 
ever, only comparison of the simulated and histor¬ 
ical daily mortality put enough constraints to the 
free data in the model (date of arrival of infected 
people and mosquito population) to allow for a se¬ 
lection of possible combinations of their values. 

As successful as the model appears to be, it is 
completely unable to produce the total mortality 
in the city, or the spatial extension of the full epi¬ 
demic. The simulations produce with BSx4 less 
than 4500 deaths, while in the historic record, the 
total mortality in the city is above 13000 cases. 
The historical account, and the recorded data, show 
that after the initial San Telmo focus has devel¬ 
oped, a second focus in the police district 13 (see 
Fig. 12) developed, shortly several other foci de¬ 
veloped that could not be tracked [5]. Unless the 
spreading of the illness by infected humans is intro¬ 
duced (or some other method to make long jumps 
by the illness), such events cannot be described. 
It is worth noticing that the mobility patterns in 
1871 are expected to be drastically different from 
present patterns, and as such, the application of 
models with human mobility |43| is not straightfor¬ 
ward and requires a historical study. 

One of the most important conclusions of this 
work is that the logical consistency of mathemati¬ 
cal modeling puts a limit to ad-hoc hypotheses, so 
often used in a-posteriori explanations, as it forces 
to accept not just the desired consequence of the 


hypotheses, but all other consequences as well. 

Last, eco-epidemiological models are adjusted 
to vector populations pre-existing the actual epi¬ 
demics and can therefore be used in prevention to 
determine epidemic risk and monitor eradication 
campaigns. In the present work, the tuning was 
performed in epidemic data only because it is ac¬ 
tually impossible to know the environmental con¬ 
ditions more than one hundred years ago. Yet, our 
wild initial guess for the density of breeding sites 
resulted sufficiently close to allow further tuning. 
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A Appendix 

i. Populations and events of the stochastic 
transmission model 

We consider a two dimensional space as a mesh 
of squared patches where the dynamics of vec¬ 
tors, hosts and the disease take place. Only adult 
mosquitoes. Flyers, can fly from one patch to a 
next one according to a diffusion-like process. The 
coordinates of a patch are given by two indices, i 
and j, corresponding to the row and column in the 
mesh. If Xk is a subpopulation in the stage k, then 
Xkihj) is the Xk subpopulation in the patch of 
coordinates (i,j). 

Population of both hosts (Humans) and vectors 
{Aedes aegypti) were divided into subpopulations 
representing disease status: SEI for the vectors and 
SEIrRTD for the human population. 

Ten different subpopulations for the mosquito 
were taken into account, three immature subpopu¬ 
lations: eggs larvae L(iy) and pupae P{i,j), 

and seven adult subpopulations: non parous adults 
Al(jj), susceptible flyers exposed flyers 

j), infectious flyers and parous adults 
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in the three disease status: susceptible A 2 s(j_j), ex¬ 
posed A 2 e(i and infectious 

The Aly y is always susceptible, after a blood 
meal it becomes a flyer, susceptible Fsyj'j or ex¬ 
posed Feyjy depending on the disease status of 
the host. If the host is infectious, becomes 

an exposed flyer Feyj) but if the host is not infec¬ 
tious, then the becomes a susceptible flyer 

Fsyjy The transmission of the virus depends not 
only on the contact between vector and host but 
also on the transmission probability of the virus. 
In this case, we have two transmission probabili¬ 
ties: the transmission probability from host to vec¬ 
tor ahv and the transmission probability from vec¬ 
tor to host avh. 

Human population Nhyy'^ was split into seven 
different subpopulations according to the disease 
status: susceptible humans Hsyyy exposed hu¬ 
mans Fley jy infectious humans FHyjy humans in 
remission state Flry jy toxic humans Flt^jy re¬ 
moved humans FIRy j-^ and dead humans because 
of the disease Hdyjy 

The evolution of the seventeen subpopulations is 
affected by events that occur at rates that depend 
on subpopulation values and some of them also on 
temperature, which is a function of time since it 
changes over the course of the year seasonally [5115] . 

ii. Events related to immature stages 

Table S) summarizes the events and rates related to 
immature stages of the mosquito during their first 
gonotrophic cycle. The construction of the transi¬ 
tion rates and the election of model parameters re¬ 
lated to the mosquito biology such as: me mortality 
of eggs, elr hatching rate, ml mortality of larvae, a 
density-dependent mortality of larvae, Ipr pupation 
rate, mp: mortality of pupae, par pupae into adults 
development coefficient and the ef emergence fac¬ 
tor were described in detail previously IHiig. 

The natural regulation of Aedes aegypti popula¬ 
tions is due to intra-specific competition for food 
and other resources in the larval stage. This regu¬ 
lation was incorporated into the model as a density- 
dependent transition probability which introduces 
the necessary nonlinearities that prevent a Malthu¬ 
sian growth of the population. This effect was in¬ 
corporated as a nonlinear correction to the temper¬ 
ature dependent larval mortality. 

Then, larval mortality can be written as: 


mlLy y +aLyj-^ X (T(ij) — 1) where the value of a 
can be further decomposed as a = with 

ao being associated with the carrying capacity of 
one (standardised) breeding site and being 

the density of breeding sites in the {i,j) patch [5]I5] . 

iii. Events related to the adult stage 

Aedes aegypti females (B1 and A2) require blood 
to complete their gonotrophic cycles. In this pro¬ 
cess, the female may ingest viruses with the blood 
meal from an infectious human during the human 
Viremic Period VP. The viruses develop within 
the mosquito during the Extrinsic Incubation Pe¬ 
riod EIP and then are reinjected into the blood 
stream of a new susceptible human with the saliva 
of the mosquito in later blood meals. The virus 
in the exposed human develops during the Intrin¬ 
sic incubation Period IIP and then begin to cir¬ 
culate in the blood stream (Viremic Period), the 
human becoming infectious. The flow from sus¬ 
ceptible to exposed subpopulations (in the vector 
and the host) depends not only on the contact be¬ 
tween vector and host but also on the transmission 
probability of the virus. In our case, there are two 
transmission probabilities: the transmission proba¬ 
bility from host to vector ahv and the transmission 
probability from vector to host avh. 

The events related to the adult stage are shown 
in Table [5] to m Table [5| summarizes the events and 
rates related to adults during their first gonotrophic 
cycle and related to oviposition by flyers according 
to their disease status. 

Table [51 and Table [3 summarize the events and 
rates related to adult 2 gonotrophic cycles, exposed 
Adults 2 and exposed flyers becoming infectious 
and human contagion. 

Table [5| summarizes the events and rates related 
to non parous adult (Adult 2) and Flyer death. 

iv. Events related to flyer dispersal 

Some experimental results and observational stud¬ 
ies show that the Aedes aegypti dispersal is driven 
by the availability of oviposition sites |44H46j . Ac¬ 
cording to these observations, we considered that 
only the Flyers Fyj^ can fly from patch to patch 
in search of oviposition sites. The implementation 
of flyer dispersal has been described elsewhere j^. 
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Event 

Effect 


Transition rate 

Egg death 



me X 

Egg hatching 

i(i.j) -t 1 


elr X 

Larval death 



7nl X -|- o X X 1) 

Pupation 

1 


Ipr X L(ij) 

Pupal death 

^(ij) ~ ^ 


(mp + par X (1 - (e//2))) X j) 

Adult emergence 

^{i,3) ~ 

1 

par X (e//2) X P(ij) 


Table 4: Event type, effects on the populations and transition rates for the developmental model. The 
coefficients are me: mortality of eggs; elr: hatching rate; ml: mortality of larvae; a: density-dependent 
mortality of larvae; Ipr: pupation rate; mp: mortality of pupae; par: pupae into adults development 
coefficient; e/: emergence factor. The values of the coefficients are available in subsections Ivi.l and ivii.l 


Event 

Effect 


Transition rate 


Adults 1 Death 



ma X Al(i j) 


I Gonotrophic cycle 



cyclel X x 


with virus exposure 

^^(i,3) 

-t + 1 

ahv 


I Gonotrophic cycle 



cyclel X Al(i 

X ({((Nh^ij-) - 

without virus exposure 


-t Fs{i,3) + 1 


-h (1 — ahv) X 

Oviposition of suscep¬ 
tible flyers 

^(ij) 

T2S(jj) 

^ ^(i,3) ~ ^ 
—3 A2S(jj) + 1 

ovr^ij) X 


Oviposition of exposed 
flyers 

^(ij) 

— t + 1 

ovr^ij) X Fe(^ijy 


Oviposition of infected 
flyers 

^(ij) 

®(ij) + 

+ 1 

ovr^ij) X Fpij^ 



Table 5: Event type, effects on the subpopulations and transition rates for the developmental model. The 
coefficients are ma: mortality of adults; cyclel: gonotrophic cycle coefficient (number of daily cycles) 
for adult females in stages Air, ahv: transmission probability from host to vector; ovri^j^ jy. oviposition 
rate by flyers in the (i,j) patch; egn: average number of eggs laid in an oviposition. The values of the 
coefficients are available in Table [T1 subsections Ivi.l Ivii.l Iviii.l and lix.l 


The general rate of the dispersal event is given 
by: P X Fyjy where /3 is the dispersal coefficient 
and Fy y is the Flyer population which can be sus¬ 
ceptible Fsyjy exposed Feyj) or infectious Fiyj) 
depending on the disease status. 

The dispersal coefficient /3 can be written as 

{ 0 if the patches are disjoint 

dif f /dfj if the patches have (3) 

at least a common point 


where dij is the distance between the centres of 
the patches and diff is a diffusion-like coefficient 
so that dispersal is compatible with a diffusion-like 
process |5]. 

V. Events related to human population 

Human contagion has been already described in Ta¬ 
ble [3 Table [3 summarizes the events and rates in 
which humans are involved. The human popula- 
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Event 

Effect 


Transition rate 


II Gonotrophic cycle 


0 - 

1 cycle2 x A2s(i x 

/A/l() X 

of susceptible Adults 2 


-t + 1 

ahv 


with virus exposure 





II Gonotrophic cycle 


0 -t - 

1 cycle2 x A2s(^ j) 

X ((((A^^qq) - 

of susceptible Adults 2 


-t + 1 


+ (1 — ahv) X 

without virus exposure 





II Gonotrophic cycle of 

A2e(i,, 

/) - 

1 cycle2 x A2e(.^ 


exposed Adults 2 


-t Fe(ij) + 1 




Table 6: Event type, effects on the subpopulations and transition rates for the developmental model. 
The coefficients are cycle2: gonotrophic cycle coefficient (number of daily cycles) for adult females in 
stages A2.\ ahv: transmission probability from host to vector. The values of the coefficients are available 
in Table [TJ subsections Ivi.l Ivii.l Iviii.l and lix.l 


Event 

Effect 


Transition rate 

Exposed Adults 2 be- 

T2e(i , 

i) 1 

[IREIP - (l/oi;r(ij))))T2e(i,^) 

coming infectious 


0 —t + 1 

Exposed flyers becom- 



{IREIP - {l/ovr(i^j))))Fe(ij) 

ing infectious 


-t A’i(ij) -1- 1 

II Gonotrophic cycle 



cycle2 x A2i^j x 

of infectious Adults 2 


-t Fif^ij) + 1 

X avh 

without human conta- 


) ~ ^ 

gion 


) -t -ffepq) + 1 


II Gonotrophic cycle 



cycle2 X 2l2j(ij) X ((((W/ipj) - 

of infectious Adults 2 
without human conta¬ 
gion 


-t Fi(ij) + 1 

+ (1 - 


Table 7: Event type, effects on the subpopulations and transition rates for the developmental model. 
The coefficients are cycle2: gonotrophic cycle coefficient (number of daily cycles) for adult females in 
stages A2; ovr(^^ jy. oviposition rate by flyers in the (f,j) patch; avh: transmission probability from 
vector to host; EIP: extrinsic incubation period. The values of the coefficients are available in Table [U 
subsections Ivi.l Ivii.l Iviii.l and lix.l 


tion was fluctuating but balanced, meaning that 
the birth coefficient was considered equal to the 
mortality coefficient mh. 

vi. Developmental Rate coefficients 

The developmental rates that correspond to egg 
hatching, pupation, adult emergence and the 
gonotrophic cycles were evaluated using the results 
of the thermodynamic model developed by Sharp 
and DeMichele m and simplified by Schoofield et 
al. gHI. According to this model, the maturation 
process is controlled by one enzyme which is ac¬ 
tive in a given temperature range and is deacti¬ 


vated only at high temperatures. The development 
is stochastic in nature and is controlled by a Poisson 
process with rate Rd(T). In general terms, Ro{T) 
takes the form 

Rz?(T) =i?D(298 K) (4) 

(T/298 K) exp((Ai7A/R)(l/298 K- 1/T)) 

"" 1 + exp(Ai7ff/R)(l/Ti/2 - 1/T)) 

where T is the absolute temperature, AH a and 
AHh are thermodynamics enthalpies characteris¬ 
tic of the organism, R is the universal gas constant, 
and Tx /2 is the temperature when half of the en¬ 
zyme is deactivated because of high temperature. 
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Event 


Effect 

Transition rate 

Susceptible flyer Death 



ma X 

Exposed flyer Death 



ma X 

Infectious flyer Death 



ma X 

Susceptible Adult 2 Death 


—>■ — 1 

ma X A2s^j 

Exposed Adult 2 Death 

^2e(j j) 

—t — 1 

ma X A2e(^ 

Infectious Adult 2 Death 


— 1 

ma X A22(j 


Table 8: Event type, effects on the subpopulations and transition rates for the developmental model. 
The coefficients are ma: adult mortality. The values of the coefficients are available in subsection I vii. I 


Event 

Effect 



Transition rate 

Born of susceptible humans 




mh X 

Death of susceptible humans 




mh X 

Death of exposed humans 




mh X 

Transition from exposed to vi- 
raemic 




(1///P) X 

Death of Infectious humans 




mh X 

Transition from infectious humans 
to humans in remission state 

Death of humans in remission state 

+ 1 

Hr(ij) - 1 



il/VP) X 

mh X 

Transition from humans in remis¬ 
sion to toxic humans 

- 1 

iyt(ij) + 1 



((1 - rar)/rP) X Pr(ij) 

Recovery of humans in remission 

Death of removed humans 

^ -f^RiP) - 1 
HR(ij) + 1 

HR(i.i) -t RR(i.i) - 1 

RR{i,j) 


(rar/rP) x 

mh X 

Death of toxic humans 



-)■ 

(mt/tP) X 

Recovery of toxic humans 

RR(i,i) +1 

RR(i,j) 


((1 - mt)/tP) X 


Table 9: Event type, effects on the subpopulations and transition rates for the developmental model. 
The coefficients are mh: human mortality coefficient; VP: human viremic period; mh: human mortality 
coefficient; IIP: intrinsic incubation period; rP: remission period; tP: toxic period; rar: recovery after 
remission probability; mt: mortality probability for toxic patients. The values of the coefficients are 
available in Table [TJ 


Table [TU] presents the values of the different coef¬ 
ficients involved in the events: egg hatching, pupa¬ 
tion, adult emergence and gonotrophic cycles. The 
values are taken from Ref. m and are discussed 
in Ref. |8]. 

vii. Mortality coefficients 

Egg mortality. The mortality coefficient of eggs is 
me = 0.01 1/day, independent of temperature in 
the range 278 K < T < 303 K |49| . 


Larval mortality. The value of oq (associated to 
the carrying capacity of a single breeding site) is 
do = 1.5, and was assigned by fitting the model 
to observed values of immatures in the cemeteries 
of Buenos Aires |8j. The temperature dependent 
larval death coefficient is approximated by ml = 
0.0H-0.9725exp(-(T-278)/2.7035) and it is valid 
in the range 278 K < T < 303 K [501152) . 

Pupal mortality. The intrinsic mortality of a 
pupa has been considered as mp = 0.01 1- 
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Develop. Cycle (0]) 

Rd{T) 

i?i>(298 K) 

XHa 

AHh 

Ti/2 

Egg hatching 

elr 

0.24 

10798 

100000 

14184 

Larval develop. 

Ipr 

0.2088 

26018 

55990 

304.6 

Pupal Develop. 

par 

0.384 

14931 

-472379 

148 

Gonotrophic c. (Al) 

cycle 1 

0.216 

15725 

1756481 

447.2 

Gonotrophic c. (A2) 

cycle 2 

0.372 

15725 

1756481 

447.2 


Table 10: Coefficients for the enzymatic model of maturation [Eq. Q]. Rd is measured in day 
enthalpies are measured in (cal / mol) and the temperature T is measured in absolute (Kelvin) degrees. 


0.9725ea;p(-(T-278)/2.7035) (MSg. Besides the 
daily mortality in the pupal stage, there is an ad¬ 
ditional mortality contribution associated to the 
emergence of the adults. We considered a mortality 
of 17% of the pupae at this event, which is added to 
the mortality rate of pupae. Hence, the emergence 
factor is e/ = 0.83 |53| . 

Adult mortality. Adult mortality coefficient is 
ma = 0.091/day and it is considered independent 
of temperature in the range 278 K < T < 303 K 

nnnii^. 


viii. Fecundity and oviposition coefficient 

Females lay a number of eggs that is roughly 
proportional to their body weight (46.5 eggs/mg) 
[551 Eg. Considering that the mean weight of a 
three-day-old female is 1.35 mg [2], we estimated 
the average number of eggs laid in one oviposition 
as egn = 63. 

The oviposition coefficient ot)r(j j) depends on 
breeding site density and it is defined as: 


opr(,j) 


6/tdep if BS(ij) < 150 
1/tdep if BS(ij) > 150 


where 9 was chosen as 0 = BS(i j)/150, a linear 
function of the density of breeding sites [3] ■ 


ix. Dispersal coefficient 

We chose a diffusion-like coefficient of diff = 830 
m^/day which corresponds to a short dispersal, ap¬ 
proximately a mean dispersal of 30 m in one day, 
in agreement with short dispersal experiments and 
field studies analyzed in detail in our previous ar¬ 
ticle jg. 


X. Mathematical description of the 
stochastic model 

The evolution of the subpopulations is modeled by 
a state dependent Poisson process FT1I57) where the 
probability of the state: 


{E, L, P,A1, A2s, A2e, A2i, Fs, Fe, Fi, 

Hs, He, Hi, Hr, Ht, HR, Hd)^i^j) 


evolves in time following a Kolmogorov forward 
equation that can be constructed directly from the 
information collected in Tables 0] to E] and in Eq. El 

xi. Deterministic rates approximation for 
the density-dependent Markov process 

Let X be an integer vector having as entries 
the populations under consideration, and ea,a = 
1... K the events at which the populations change 
by a fixed amount Aq, in a Poisson process with 
density-dependent rates. Then, a theorem by 
Kurtz allows us to rewrite the stochastic pro¬ 
cess as: 


X{t) = A(0) -f 

Ct^l 


f 

0 




( 6 ) 


where a;c[(A(s) is the transition rate associated 
with the event a and Y(x) is a random Poisson 
process of rate x. 

The deterministic rates approximation to the 
stochastic process represented by Eq. m consists 
of the introduction of a deterministic approxima¬ 
tion for the arguments of the Poisson variables Y{x) 
in Eq. (E]) |M1[5H] . The reasons for such a proposal 
is that the transition rates change at a slower rate 
than the populations. The number of each kind of 
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event is then approximated as independent Poisson 
processes with deterministic arguments satisfying a 
differential equation. 

The probability of ria events of type a having oc¬ 
curred after a time dt is approximated by a Poisson 
distribution with parameter Aq,. Hence, the proba¬ 
bility of the population taking the value 

K. 

X = Xo + Y, (7) 

a=l 

at a time interval dt after being in the state Xq is 
approximated by a product of independent Poisson 
distributions of the form 


Probability(m ... n^, dt/Xg) = Pa{K) (8) 

a—1 

and 

\na 

P:,...nM=eM-^c.)^ (9) 

Tloi'. 

whenever X = Xq + X]a=i has no negative 

entries and 

°° y 

^“i...nE(Aa) = exp(-A„) 

i—Ua 

rio-l y 

= 1 - exp(-Ac) Y ^ (10) 

if {rii} makes a component in X zero (see Ref. m) 
Finally, 

dXa/dt =< LUaiX) > (11) 

where the averages are taken self-consistently with 
the proposed distribution (Aa(0) = 0). 

The use of the Poisson approximation represents 
a substantial saving of computer time compared 
to direct (Monte Carlo) implementations of the 
stochastic process. 
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